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PACS 64 . 70 . k j - Glasses 
PACS 64 . 70 . pm - Liquids 
PACS 64 . 70 . pv - Colloids 

Abstract. - We show that the dynamics of soft-sphere systems with purely repulsive interac- 
tions can be described by introducing an effective hard-sphere diameter determined using the 
Andersen- Weeks-Chandler approximation. We find that this approximation, known to describe 
static properties of liquids, also gives a good description of a dynamical quantity, the relaxation 
time, even in the vicinity of the glass transition. 



Relaxation in a liquid proceeds by numerous pathways 
that each depend on temperature and pressure. It has 
remained a challenge to delineate different modes of re- 
laxation and to bring them into a framework where they 
can be treated on a common footing [IHl] . To reach such 
a goal, it is necessary to understand what combinations 
of variables control liquid dynamics in different regimes of 
temperature and pressure. 

Some insight was gained by expressing the liquid relax- 
ation time, T, in units of a time scale, \fmjpu, defined by 
the particle mass, ra, diameter, a and pressure, p. We note 
that this is not the most obvious choice for expressing the 
time scale; the conventional choice is ^mtr^/e, where e is 
the energy scale of the inter-particle potential. The choice 
yjm/pcr focuses on the elementary time for a particle re- 
arrangement event to occur driven by the pressure. The 
dimensionless relaxation time can be written as a function 
of dimensionless variables [TOl : 



Tyjpa/m = f{T/pa^,pa^/e, ..), 



(1) 



where T is the temperature (or thermal energy in units 
where Boltzmann's constant is set to unity). Other pos- 
sible dimensionless variables include the polydispersity, 
or quantities describing the shape of the particles or the 
potential. It was shown that if we confine ourselves to 



monodisperse spheres with finite-ranged repulsions then 
there is a major simplification [TU] when per' je is small 
so that there is little particle overlap. In that regime, 
the relaxation times for many systems with different in- 
teractions can be collapsed onto a single master curve that 
corresponds to the hard sphere limit [10 . Thus, 



lim 



-\fpajm = fnsiT/pcr^). 



(2) 
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In this regime, the simple interpretation can be made that 
relaxation in the hard-sphere limit proceeds by thermal 
energy working against pressure to open up free volume. 

Unlike hard spheres, soft spheres can rearrange by cross- 
ing the potential barriers imposed by their neighbors. 
Once soft spheres begin to overlap, the dynamics deviate 
from those expected in the hard-sphere limit and depend 
sensitively on the interaction potential, as one might ex- 
pect [10) . Here we analyze the situation where the particle 
overlap is no longer negligible so that soft particles do not 
behave simply as hard spheres with the same diameter, 
as in eq. [H Instead, the system responds as if it were 
made up of particles with a smaller effective hard-sphere 
diameter, Coff- 

It has long been known that the static properties of 
soft repulsive spheres can be very well approximated by 
hard spheres with an appropriately chosen effective di- 
ameter [TTHT5] . For damped soft-sphere systems, self- 
consistent generalized Langevin theory has been used to 
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generalize these static results to the dynamics [7]. 

Here, we calculate the effective diameter using the 
procedure proposed by Andersen, Weeks and Chan- 
dler (AWC) [TS] for the static properties, as Guevara- 
Rodriguez, et al. have done for damped spheres. We 
replace the relaxation time, r of a liquid of soft spheres 
at temperature, T, and number density, p, with an effec- 
tive relaxation time Tcs of a liquid of hard spheres with 
effective diameter, cTcff, and pressure, Peff. This hard- 
sphere system is constrained to have the same values of 
p, T and to. The hard-sphere equation of state relates 
the packing fraction, 0, to T/pa^. Thus, the packing frac- 
tion of the effective hard-sphere system, t/ieff = 7rpCTJ^ff/6, 
corresponds to a certain value of T/pescr^ff- This speci- 
fies the relaxation time of the effective hard-sphere liquid: 
Teffx/pcff gcff/m = /Hs(r '/PcffO-j^ff)- Our central finding is 
that Ty/pa/rn = cTcffv/PcffO'off/'Ti, where c = 1 in case 
of finite-range repulsions. For long-ranged interactions, 
where the particle diameter is not well-defined, c depends 
on the potential but is independent of temperature and 
pressure. Therefore the relaxation time of the soft-sphere 

system is 

T = CTcff v^PcffCcff /per (3) 

where Peff and Ueff are the pressure and diameter of the 
hard-sphere system given by the AWC approximation and 
p and (J are the pressure and diameter of the original soft- 
sphere system. We find that Eq. [3] describes the relaxation 
time reasonably well for many different repulsive poten- 
tials over a wide range of temperatures and densities, in- 
cluding near the glass transition. 

System and simulation methods. We consider 
a three-dimensional system of monodisperse spheres that 
interact according to a pair interaction Vijij) where rij is 



the distance between particles i and j: 



V{nj) 




(7 } 
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with different exponents a > 0. We also consider the 
Weeks-Chandler- Anderson potential [16] 
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and power-law repulsions. 
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(5) 



(6) 



with /3 = 6 or 12. In the simulations, the interactions 
in Eq. [6] were cut off at 2.5cr and a correction term cor- 
responding to the force at the cutoff is subtracted from 
Vifij) such that the first derivative of the potential em- 
ployed in the simulations is continuous. 

To calculate the relaxation time, we performed molecu- 
lar dynamics simulations at fixed pressure p and temper- 
ature T of systems of iV = 1000 particles, using periodic 
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Fig. 1: Dimensionless relaxation time ry^pa/m, for different 
potentials as defined in Eqs.[4]-WCA, as a function of pressure 
po-^/e for T/{pG^) = 0.06 (blue) and 0.2 (red). As a increases, 
the dimensionless relaxation time decreases, as shown by the 
arrows labelled a. 



boundary conditions. For the case of hard spheres, a = 
in eq. HI we used an event-based algorithm. We define 
the relaxation time r as the time where the mean square 
displacement is ("'"^(t)) — a^ . We estimated the error 
bar in the relaxation time by averaging over 20 different 
starting conditions. The standard deviation among these 
runs divided by 20 yields an error bar of about 1.4% of 
the relaxation time for several different state points for 
the different potentials studied. At each temperature and 
pressure, we also calculated the mean number density, p. 
In order to detect any crystallization of our samples, 
we calculated the bond-orientational order parameter Qe 
[17]. For our soft-sphere systems, Qg never exceeded 0.27 
which is well below the value of a face-centered-cubic crys- 
(4) tal (Ql^'^ « 0.575 [H]). The magnitude of Qg is compa- 
rable to values measured for the large spheres in bidis- 
perse systems at similar relaxation times and pressures. 
We observed no increase of Qq with simulation time. In 
addition, the third-order invariant has values We ~ 0.05, 
very close to the values reported for Lennard- Jones super- 
cooled liquids and considerably higher than for crystalline 
clusters [T7]. These results indicate that there is no ap- 
preciable crystallization in our samples. 

Analysis and Results. — Figure [1] shows the dimen- 
sionless relaxation time versus pa^ j e for several different 
potentials at two values of T/pa^ . For fixed T/pa^, the 
relaxation time approaches its common hard-sphere value 
(shown in fig. [1] as gray horizontal lines) at small per" je 
for all potentials shown, in agreement with previous re- 
sults [To]. However, as pcr"^/e increases, the relaxation time 
decreases in a potential-dependent fashion. The results 
shown extend to pressures of order O.le/cr''; for typical 
Lennard- Jones parameters for nonpolar molecules, where 
ejU - 100 - imK and cr ~ 3 - 61 [151, this corresponds 
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to pressures on the order of several tens of atmospheres. 

The first step in understanding the decrease of relax- 
ation time in fig. [1] with pressure is to define an effec- 
tive hard-sphere diameter to describe the liquid. Rowlin- 
son [mill] proposed to define a^s based on an expansion 
of the partition function in powers of an inverse steep- 
ness of the potential. This approach was further devel- 
oped by Barker and Henderson [13] and Andersen, Weeks 
and Chandler (AWC) [15]. The idea of the AWC approx- 
imation is to consider a set of potentials that interpolate 
smoothly between the effective hard-sphere potential with 
diameter tToff and the original potential V{r). The free 
energy of the original system is then expanded around the 
hard-sphere limit in a functional Taylor expansion. AWC 
chose the effective hard-sphere diameter a^e such that the 
first order term in the Taylor expansion vanishes. This 
yields the condition 



drycffir) 



exp(-F(r)/r) -exp(-V^eff(r-)/T) 



-0 (7) 



where T4ff (r) is the potential for hard spheres of diameter 
tJoff and ycsir) is the hard-sphere cavity pair correlation 
function defined by j/cff(?') = ffoff('') exp(T4ff(r)/T), where 
gesif) is the pair correlation function of the hard-sphere 
system. 

We solve eq. [7]to obtain the AWC effective hard-sphere 
diameter, aeS- For ycs{r), we employ the polynomial ap- 
proximant introduced by Grundke and Henderson [19j . 
which uses the Verlet-Weis-approximation [20] for r > a^e, 
where j/cff(^) = 5cff, and extends it to r < cFcS- The 
Grundke- Henderson approximation to j/eff has been proven 
to be in excellent agreement with simulation results for 
hard spheres [21]. 

The ratio of the resulting effective diameter to the origi- 
nal soft-sphere diameter, fToff/c, is plotted versus pa^ /e in 
fig. [2] for several different potentials. At fixed T/pa^, the 
ratio decreases with increasing pa^ /e. This is expected, 
since the overlap increases with increasing pa^ /e. Addi- 
tionally, fig. [2] shows that the effective hard-sphere diam- 
eter decreases with increasing T/pa^ at fixed pa^ /e. This 
is also expected, since higher overlaps, which cost more 
energy, are more probable at higher temperature. 

While the effective hard-sphere system and the original 
soft-sphere system have the same temperature and num- 
ber of particles, they have different pressures. From the 
effective packing fraction of the hard-sphere system, 0off, 
we calculate Pesc^g/T, using the Carnahan-Starling ap- 
proximation [35]. Thus, for each value of T and p studied 
in the original soft-sphere system, we extract a packing 
fraction (f)cs and T/pcscr^g for the effective hard-sphere 
system. 

For finite-ranged repulsions as in Eqs.|4][5l the potential 
unambiguously defines the particle diameter a. However, 
for power- law repulsions, Eq.lHJ the choice of a is arbitrary. 
To compare our results for these different potentials, we 
introduce the parameter c. For finite-ranged repulsions we 
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Fig. 2: The ratio of the effective hard-sphere diameter to the 
soft-sphere diameter, (Jcff/cr, calculated from the Andersen- 
Weeks-Chandler (AWC) approximation 15 as a function of 
pa^/e at T/pa^ = 0.06 (blue, solid lines) and T/pa^ = 0.2 
(red, dashed lines). Lines are added to connect data points 
for the same potential at the same T/pa^. Different symbols 
correspond to different potentials as defined in the legend of 

fig.m 



have c = 1 while for power-law interactions we determine 
c by a fit to Eq. |3] that depends on the potential but does 
not depend on temperature or pressure. We find c = 1.31 
for /3 = 6 and c = 1.85 for l3 = 12. 

The relaxation times for the different potentials are 
shown as a function of T/pa^ in fig. EJ^a). Fig. ^ih) shows 
the dimensionless relaxation time T^pa/m/c versus effec- 
tive packing fraction 4>es'^^ (inset) and versus T/pcS<J^Q 
(main figure). These figures show that the data collapse 
fairly well onto the hard-sphere curve /hs (gray line) of 
eq. [U Thus, /hs gives a reasonably good description of 
the relaxation time of all the systems with repulsions that 
we have studied. This is a remarkable result. It shows 
that for a given soft-sphere system, one can calculate an 
effective ratio of temperature over pressure to obtain a 
reasonable estimate of the dynamics, even close to the dy- 
namical glass transition. 

The collapse shown in fig. Eljb) is best where the ratio 
Ocs/o is sufficiently close to unity. Large symbols in fig. [3] 
denote data where CTcs/ct > 0.9 while small symbols are 
used for systems with higher pressure where Ueff/cr < 0.9. 
Some of the data points with small Ucs/ a are systemati- 
cally above the hard-sphere curve (especially at high tem- 
peratures, T/pesa^ff > 0.2), i.e., the relaxation time is 
systematically larger than that of the corresponding hard- 
sphere system. For such small CTcff/c values, overlaps of 
three or more spheres become significant, leading to cor- 
rections to the AWC approximation and a breakdown of 
collapse. Note from Fig. [2] that all of the data for the 
WCA potential fall in the range where fToff/c > 0.9 and 
are represented by large symbols; it is only for unphys- 
ically soft repulsions that there is substantial overlap at 
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Fig. 3: (a) The dimensionless relaxation time T^/paJrn versus 
T/pa^ for soft-spheres for different potentials and values of 
■pcP' jt. (b) Dimensionless relaxation time T^JpaJm/c versus 
the effective volume fraction (/)cff (inset) and versus T/poffCToff 
(main figure). For finite-ranged interactions the symbols for 
the data points are the same as in the legend to fig.[T] Different 
colors correspond to different values of T/pa'^: T/pa^ — 0.04 
(black), 0.05 (purple), 0.06 (blue), 0.08 (cyan), O.fO (green), 
0.15 (orange) and 0.20 (red). Data points where (Jas/cr > 0.9 
are shown with large symbols, data with smaller (Toff/cr are 
represented by small symbols. Arrows indicate the direction of 
increasing pressure at each value of T jprj''' . Grey data points 
were obtained for other values of T jpa^ and are only shown 
for (Jes/<y > 0.9. Solid triangles and circles denote data with 
power-law interactions, using /? = 6 and /3 = 12, respectively. 
The typical statistical error is denoted in (b) by two yellow 
points. 



high pressures so that (Tcff/cr < 0.9. 

Because it is easier to obtain data points with high 
values of T\Jpa/m with simulations of soft spheres at 
higher pressures than for the corresponding hard-sphere 
system, we can use soft-sphere calculations to extend the 
relaxation-time curve for hard spheres. This is shown 
in fig. [21 where the gray curve, representing our hard- 
sphere simulations, ends at T^ppajm « 10^ due to the on- 
set of crystallization but soft-sphere results extend above 



T\J-pa jm — 10^. Note, however, that the extension of 
the master curve /hs to higher relaxation times relies on 
knowledge of the metastable liquid branch of the hard- 
sphere equation of state. In particular, we used the 
Carnahan-Starling approximation to obtain fig. [3l^b); this 
approximation must fail at sufficiently high pressures. 

Discussion. — We note that the scaling collapse 
shown in Fig.[3]is different from the one found by Berthier 
and Witten [UIl] , who studied bidisperse mixtures of soft 
spheres. In their analysis, the data are collapsed onto two 
scaling curves, one for 4> < (po and the other for > (/)q, 
where 0o and the two scaling exponents are chosen to ob- 
tain the best collapse. By contrast, the scaling collapse 
based on the AWC approximation collapses all of the data 
onto a single curve, rather than two curves. The data 
shown in fig. [3] span a much broader range of packing frac- 
tions and temperatures than the data included by Berthier 
and Witten. In addition, we have included data for many 
different potentials, whereas Berthier and Witten's analy- 
sis was only for the case a = 2 in eq. HI For finite-ranged 
interactions, our analysis yields a reasonable collapse with 
no adjustable parameters. For long-ranged interactions, 
there is only one free parameter, c, that depends on the 
arbitrary choice of the original diameter a. 

The AWC approximation leading to the collapse shown 
in fig. [3] also differs from collapses of the dynamics based 
on the excess entropy. Based on the ideas of Rosenfeld 
[9] and Gnan et al. [3j , a number of studies have explored 
the connection between the dynamics and the excess en- 
tropy [IHS]. The ratio T/pcfta^g could perhaps be viewed 
as a proxy for the excess entropy; for hard spheres the 
excess entropy depends on T/pcsa^g. At least for some 
potentials [6], however, what appears to collapse with ex- 
cess entropy, is the diffusivity made dimensionless by the 
energy scale of the potential, e. This is fundamentally dif- 
ferent from the relaxation time made dimensionless by the 
pressure; the difference involves the ratio pcr^/e, which is 
not fixed. 

In summary, we have extended the results of Xu, et 
al- [To] by showing that the relaxation times of liquids 
composed of spheres interacting via finite-ranged repul- 
sions can be mapped up to moderately high pressures 
onto a system of hard spheres with an appropriate effec- 
tive diameter: Ty^pa/m — cfusiT/pa^) where c = 1 for 
finite-ranged repulsions. This is accomplished by replac- 
ing the soft-sphere system with a system of hard spheres 
of smaller diameter, a^g, where Ueff is the same value 
that provides an excellent approximation to structural and 
other static properties |15j . This "principle of dynamic 
equivalence" was first noted by Guevara-Rodriguez and 
Medina-Noyola [7l[8] in the context of Brownian particles, 
and has been generalized in this paper to undamped sys- 
tems with a wide range of repulsive potentials. 

The relaxation time for soft spheres is fairly well de- 
scribed by the hard-sphere curve up to moderate pres- 
sures. This suggests that the physics that controls the 
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glass transition for hard spheres also controls it for soft- 
sphere systems. This is a striking affirmation of the util- 
ity of the free-volume concept for soft-sphere liquids. The 
free- volume approach to dynamics in liquids, which has 
been applied to fluids ranging from polymers [531[2S] to 
aqueous electrolytes IJB] to glassforming liquids [37], as- 
sumes that holes must open up in a dense liquid in order 
for molecules to rearrange. In liquids of hard spheres, the 
concept of a hole is well-defined and the free- volume distri- 
bution can be calculated explicitly and efficiently [28]. In 
a system of soft spheres, however, the definition of a hole, 
and consequently of free volume, has been unclear. The 
replacement of soft spheres with hard spheres of smaller 
diameter, as we have done here, resolves this uncertainty. 
Our results suggest that free volume is controlled by two 
distinct mechanisms in liquids of repulsive spheres. First, 
thermal fluctuations perform work against the pressure to 
open up free volume. This mechanism is governed by the 
quantity T/pa^. Second, at a given T/pa^, the pressure 
forces soft spheres to overlap so that they behave like hard 
spheres of smaller diameter. This mechanism is controlled 
by pu^ /t. These two mechanisms combine to determine 
the ratio of thermal energy over an effective pressure of 
the system. The collapse of all of our soft-sphere data 
onto a hard-sphere relaxation time curve implies that the 
relaxation time of the soft-sphere liquid depends largely 
on this ratio. The next major challenge is to test whether 
systems with attractive as well as repulsive interactions 
can be described by generalizations of the two mechanisms 
for opening up free volume identified here. Furthermore it 
would be interesting to see whether our method also works 
in systems with shear, where temperature over pressure 
has proven to be a useful control parameter, especially in 
the limit of small pressures [29] . 
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